Dynamic Regimes in Films with a Periodic Array of Antidots 
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' We have studied the dynamic response of Pb thin fihns with a square array of antidots by means 

of ac susceptibihty x{T, H) measurements. At low enough ac drive amphtudes h, vortices moving 
, inside the pinning potential give rise to a frequency- and /i-independent response together with a 

^ ■ scarce dissipation. For higher amplitudes, the average distance travelled by vortices surpasses the 

' pinning range and a critical state develops. We found that the boundary h*{H,T) between these 

regimes smoothly decreases as T increases whereas a step-like behavior is observed as a function of 
^ . field. We demonstrate that these steps in h*{H) arise from sharp changes in the pinning strength 

corresponding to different vortex configurations. For a wide set of data at several fields and temper- 
^ I atures in the critical state regime, we show that the scaling laws based on the simple Bean model 

P5 . are satisfied. 

o ■ 

■ 

5^ \ I. INTRODUCTION 

^ r AC susceptibility X = x' + measurements are a powerful tool widely used for investigating the vortex dynamics 
^ ; ' in the mixed state of Type-II superconductors. From the experimental point of view, this inductive technique has the 
, advantage of being inexpensive, highly sensitive, and very simple. On top of that, it allows one to access and probe 

■ different vortex regimes and phases by changing either the "sensing parameters", like frequency /, amplitude h and 

1 ' wave shape of the alternating excitation, or the thermodynamic variables like the external field H and temperature 
■ O T. For example, in the vortex soHd phase, small values of the external excitation h induce an oscillatory motion of 
^ \ the flux lines inside the pinning potential characterized by a very small dissipation and a response x independent of 
2 . /i.Pl In this so-called Campbell regime, it is possible to determine the average curvature of the pinning potential a 

assuming a parabolic potential well. 0, Q 

At higher amplitudes h, the vortex displacement approaches the pinning range and the pinning potential departs 
from the harmonic approximation. At this stage, a nonlinear response is observed and the screening x' becomes h- 
and /-dependent. Eventually, at high enough amplitudes, the average distance travelled by vortices may surpass the 
distance between nearby pinning centers (inter- valley motion) and a critical state (CS) develops. 
\ Unlike the Campbell regime, where the vortex configuration is preserved after measuring, in the critical state regime 
fSj ' the initial vortex configuration may be strongly altered during the measurement process. However, it is precisely this 
' invasive character of the measurements in the CS which allows one to determine the critical current density J (i.e. 
\l the depth of the potential well) , thus providing complementary information to that obtained in the Campbell regime 
(i.e. the curvature of the pinning potential). 

Typically the transition between the Campbell and the critical state regimes is broadened by topological and 
energetic disorder. These effects can be substantially reduced by introducing a regular array of pinning centers. In 
a recent workQ we have shown that this reduction of disorder leads, within the Campbell regime, to well defined 
. transitions in the ac response at the commensurability fields H = Hn = n^o/cP, where $o is the superconducting 
' fiux quantum and d is the period of the square pinning lattice. 
O 1 It is generally thought that the vortex distribution in this kind of samples consists of terraces of homogeneous 
O vortex distribution connected by narrow walls where the vortex density changes abruptly. Q Within this picture, 
the applicability of standard critical state models with smooth field profiles is highly questionable, and even more 
doubtful is the validity of the simple Bean critical state model commonly used to extract the critical current from 
. , magnetization measurements. However, the lack of any corroborative evidence of such terraced critical state profile 

■ in samples with periodic pinning stimulates further experimental studies to gain insight into this rather unexplored 
" " topic. 

In this work we study the ac dynamic response of Pb thin films with a square pinning array of holes. We determine 
the crossover field from linear to nonlinear regimes as a function of temperature and field. We found that this 
boundary smoothly decreases with increasing temperature and exhibits a step-like field dependence with clear jumps 
at every matching field. For large amplitudes, in the critical state regime, the scaling relations derived from a simple 
Bean model are accurately satisfied for a large range of fields and temperatures. This analysis allows us to estimate 
quantitatively the temperature and field dependence of the critical current density as well as to determine the onset 
of the critical state regime. 
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II. EXPERIMENTAL DETAILS 

The experiments were conducted on two Pb thin films with a square antidot array of period d ~ f.5/im, which 
corresponds to a first matching field Hi = 9.2 G. The antidots have a square shape with a size b — 0.8/im. The name, 
critical temperature and thickness S of the used samples are, AD15 {6 = 13.5 nm, = 7.1 K) and AD65 {S = 65 nm, 
Tc = 7.21 K), respectively. From the Tc{H) slope we have estimated, for both samples, a superconducting coherence 
length ^(0) - 33 ± 1 nm. 

The predefined resist-dot patterns were prepared by electron-beam lithography in a polymethyl metacrylate/methyl 
metacrylate (PMMA/MMA) resist bilayer covering the Si02 substrate. The use of two different resist layers is 
necessary to obtain an overhang profile which eventually guarantee that there will be no contact between the material 
deposited on the substrate and the resist layer. A Ge(20 A)/Pb/Ge(200 A) film was then electron-beam evaporated 
onto this mask while keeping the substrate at liquid nitrogen temperature. Finally, the resist was removed in a liftoff 
procedure in warm acetone. 

The ac-measurements were carried out in a commercial Quantum Design-PPMS device with drive field amplitudes 
h ranging from 2 mOe to 10 Oe, and the frequency / from 10 Hz to 15 kHz. In this window of frequencies we have 
found that x is weakly dependent on / and therefore we report results obtained at a single value / = 3837 Hz. In all 
cases, the data were normalized by the same factor corresponding to a total step Ax' = 1 with H —0. 

III. BOUNDARY OF THE LINEAR REGIME 

In order to identify the different dynamic regimes, we carried out a series of measurements of x = x' + *x" 
a function of the ac excitation h for several temperatures at fixed dc field H. Some of these measurements for the 
sample AD15 at H =5 Oe, are shown in figure ^ These curves clearly indicate the existence of a Campbell regime at 
low h values, characterized by a susceptibility independent of h (see dashed lines in the figure) and x" ^ 0. At higher 
amplitudes, the dissipation departs from zero and a crossover to an amplitude-dependent nonlinear regime takes place 
(see black arrows in this figure). 
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FIG. 1: Ac-susceptibility x ~ x' + *x" ^ ^ function of the drive field h for several temperatures at H = 5 G. Horizontal dashed 
lines indicate the range of the linear response and the arrows the onset of the nonlinear regime. 

We use the criterion x" — 0.01 as a reliable estimation of the onset of the nonlinear regime. Performing this analysis 
for different temperatures and dc fields we determine the boundary h* (T, H) between the Campbell and the nonlinear 
regime. The obtained h*{t — T /Tc) for H — b G, \s shown in the main panel of Figure|21as open symbols. 

For drive fields h < h* , the pinning potential can be approximated by a parabolic well and the ac-response arises 
from vortices performing small intra- valley oscillations inside their pinning site. For h ^ h*, the average vortex 
displacement approaches the pinning range where the harmonic approximation breaks down and a nonlinear response 
is detected. Eventually, at higher amplitudes, vortices are able to hop out the potential well and a critical state 
develops. 
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Since the limits of the hnear regime are given mainly by the strength of the pinning centers, it is expected that the 
extension of this regime decreases as temperature increases, in agreement with the observed behavior. Interestingly, 
in these patterned samples, we are able to substantially change the nature of the pinning potential, and thus its 
strength, by simply changing the external magnetic field. Indeed, for fields H < Hi the density of pinning centers 
is larger than the density of vortices and each flux line is strongly pinned by an antidot. For H > Hi extra vortices 
could still be attracted by the antidots if the maximum number of flux lines that an antidot can hold ris > 1. In 
this case, since multiquanta vortices are weaker pinned than a single vortex, a reduction in the extension of the linear 
regime is expected at each matching field. For fields H > rig Hi the incoming vortices will sit in interstitial positions 
caged by the repulsive interaction of their strong pinned neighbors. This family of vortices, much weaker pinned and 
with higher mobility, 's'l will dominate the ac- response and a more dramatic reduction of h* should occur. 




FIG. 2: Main panel: dynamic h — T phase diagram at = 5 G. The h*{T) curve represents the onset of the nonlinear response 
and h'^''{T) indicates the beginning of the critical state regime. Inset: dynamic diagram in the h — H plain at fixed temperature. 
The lines are guides to the eye. 

This picture is consistent with the behavior of h*{H) shown in the inset of Figure |21 at T = 7 K = 0.985). In 
this figure we observe that for < < the h* data scatter around 21 mG (indicated by a straight line in the 
figure) whereas for Hi < H < H2 the average value drops to 18 mG. In both field ranges, the data points exhibit 
a similar dispersion showing small peaks near Hi/2 and 3i?i/2, where the interstitial flux line lattice form a highly 
stable configuration. The small step in h*{H) at Hi indicates the transition from single quantum pinning to double 
quanta pinning. At H ^ H2, h* undergoes a stronger reduction down to 12 mG, followed by a different field evolution 
as H increases. This more pronounced jump in h* suggests that for H > H2, a further increase of the occupation 
number per hole is no longer energetically convenient and hence, incoming vortices will locate in interstitial sites (i.e. 

A first theoretical estimation of the maximum number rig of vortices trapped in a cylindrical cavity was done more 
than 30 years ago by Mkrtchyan and Schmidt 6i|. According to this calculation, Ug « &/4^(T). In our particular case, 
this gives rig ~ 1. Although this value is smaller than the experimentally determined above, it is important to note 
that this model underestimates the real since it considers a single vortex line interacting with a cylindrical insulating 
defect with radius b ^ X. Clearly, this hypothesis is not fulfilled in our system where b ^ 2A(r). Additionally, it 
is also expected that increases as the applied dc field is increased. 0) H 13 An extension of the original work of 
Mkrtchyan and Schmidt to arbitrary large cavity radius has been recently done by Nordborg and Vinokur^3| using 
the London approximation. Similar results were found by Buzdin'P| using the method of image vortices. This author 
showed that in a triangular vortex lattice, a two-quanta vortex becomes energetically favorable for temperatures such 
that b^ < ^(T)A(r)^, a condition that, in our sample, is satisfied a,t t < 0.995, in agreement with the experimental 
result. 



IV. CRITICAL STATE REGIME 



As we pointed out above, for h > h*, the vortex motion is no longer restricted to local oscillations around the 
pinning centers and vortex excursions may be larger than the distance between nearby pinning sites. The ac response 
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in this case can be described within a critical state model. In the simplest scenario of a Bean model, the persistent 
current density J is uniform and close to the critical value Jc, everywhere in the sample where the inter- valley motion 
takes place. In this regime, unlike the Campbell regime, the field penetration A increases linearly with h as.[ll.] 



It is important to note that the ac response x is solely determined by the penetration depth A of the ac excitation. 
As a consequence, it is possible to find, in the critical state regime, a set of data at different T and /i, such that 
produce the same penetration A, and therefore the same x- For example, according to equation if all data in 
figure n were in a well developed critical state regime, curves measured at constant T should collapse on a single 
one when, for each curve, the horizontal axis h is scaled by a suitable factor J[T)^^. The resulting curve from this 
procedure is the Bean penetration depth A(x') « h/J{T). This scaling is shown in the main panel of FigureOfor 
temperatures ranging from 6 K to 7 K. We observe a good overlap at low screening values < 0.7) and a break 
down of the scaling when approaching to the Campbell regime at higher x' values, where the field penetration and 
therefore x', become /i-independent. 




FIG. 3: Screening x' bjs a function of the Bean penetration depth A obtained by scaling the curves x'C*) showed in Figure Q 
The upper inset shows the critical current as a function of temperature obtained from different methods for a plain film and a 
patterned film at J? = 5 G. The lower inset shows the efficiency e{T) — J{antidots)/J{plainfilm). 

The set of parameters J(T) chosen to obtain the superposition of the low screening data, determine the temperature 
dependence of the persistent current as is shown by star symbols in the upper inset of figure 13 Alternatively, the 
persistent current J{T) can be also obtained from x"W curves measured at fixed T, by tracki ng t he position /i™°^(T) 
of the maximum dissipation. Indeed, exact calculations for the particular case of a disk showpd| that x" maximizes 
when A — 5 and thus from equation ^ we have that for each T, J(T) « K^"-^ /5. This is shown as solid circles in 
the upper inset of figure |21 A similar argument can be applied to obtain J{T) from the temperature at which the 
dissipation maximizes in x" {T) curves. 12] The result of this procedure is shown by open squares in the same inset. 

For comparison, in the upper inset of figure|21we also show J(T) for a Pb plain film without antidots deposited 
simultaneously with the patterned film. The efficiency s{T) = J {antidots) / J {plain film) of the antidots as correlated 
pinning centers with respect to the random intrinsic pinning of the reference film, is shown in the lower inset of 
figure 131 We observe that e monotonically decreases as T decreases from Tc- This result is in agreement with the 
experimental fact that commensurability effects originated in the periodicity of the pinning array, progressively fade 
out as intrinsic pinning becomes more relevant at lower temperatures. The very high e values observed in this sample 
can be taken as a confidence test of the quality of the grown films. 

Let us now analyze the crossover from the nonlinear intermediate regime to the critical state regime. According 
to the theory, in the simplest case of a cylinder in a parallel field, the critical state regime arises when the Campbell 
penetration depth A(r) becomes equal to the Bean penetration A{T,h), i.e. when h ~ J(r)A(T).^3| In a previous 
work*^ we have determined for a similar sample, that \{t = 0.90) ~ 250 nm and using the J{T) from figure 13 
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we obtain a crossover field /i"^** ~ 40 G, which turns out to be far above our experimental estimation of h'"^ '--^ 1 G 
according to figure 01 Recently, Pasquini et al.0 have pointed out that in the case of transverse geometry, since 
the actual range of field penetration does not coincide with the penetration depth, the previous criterion should be 
modified. Assuming that the range of field penetration in both regimes coincides when \x'\ = 0.5, the authors show 
that the crossover should occur at /i'^" ~ yj%5 /2R\{T)J{T), where R is the radius of a disk shaped sample. In our 
case, the geometric factor ^JQ6/2R ^ 0.007 shifts the crossover field down to /i^* ~ 0.3 G. Although this crossover 
field is somewhat smaller than the experimental value, we should note that whereas this simplified model correctly 
accounts for the reduction of the transition field /i"^^, it still involves a certain degree of arbitrariness since different 
criteria lead to different geometric factors. 

In fact, figure^ shows that the scaling derived from the critical state becomes valid for \x'\ < 0.7 instead of 
Ix'l — 0.5. From the penetration depth at this screening value. Across 6 nm, and the obtained critical current J{T), 
we can estimate the onset of the critical state regime /i"^* — Across J{T). This boundary is represented by solid circles 
in the main panel of Figure [3 It is apparent in this figure that the transition between the linear and the critical state 
regime is not a sharp crossover but instead, a broad intermediate non linear regime lies in between. 

There are two different sources that give rise to this crossover regime. First, it appears as a natural consequence of 
the anharmonicity of the pinning potential. Second, due to the inhomogeneous current distribution, nonlinearities first 
appear at the border of the sample and then, as h increases, the boundary separating inter- and intra-valley vortex 
motion, moves towards the center of the sample. This coexistence between linear and nonlinear regime gives rise to a 
nonlinear response. An additional effect that may further expand the intermediate regime is the presence of disorder 
in the distribution and energy of the pinning centers. However, as we pointed out above, in the studied samples owing 
to the regular array of pinning centers, the topological and energetic disorder is reduced down to minimum levels and 
should not play a significant role. 

A similar scaling analysis can be performed to obtain the field dependence of the critical current. To that end, we 
have measured the ac susceptibility response of the AD65 sample, as a function of h, for several applied dc fields H 
a.t T — 7 K (t = 0.97). These curves are shown in the lower inset of figure^ for fields ranging from G to 56 G. 
Strictly, equation is valid within the Bean critical state scenario, i.e. with a critical current Jc independent of 
g.p^ This shortcoming can be avoided either by applying an ac drive much smaller than the dc field, or provided 
that Jc is insensitive to the small field oscillations. 
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FIG. 4: Screening x' a function of the Bean penetration depth A obtained by scaling the curves showed in the lower inset. 
The solid line depicts the theoretical expectations according to the Bean model. The lower inset shows the real part of the ac 
response x' a-s a function of the amplitude h for several _ff at T = 7 K (t = 0.97). The upper inset shows the critical current 
as a function of field at t = 0.97, obtained following the procedures explained in the text. 

In the main panel of figure 01 we show the results of a scaling procedure similar to that previously explained to 
obtain A(x')- For \x'\ < 0.6, a fairly good overlap of all data points is clearly observed. This result indicates that in 
this field range, equationlQ]) is satisfied and therefore the Bean model applies. The solid line depicts the theoretical 
expectations according to the Bean model for a thin disk.^J If is worth to note that this is not a fitting curve but a 
theoretical prediction where no free parameters are involved. 
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The critical current J{H) derived from this scahng is shown in the upper inset of figure 0] together with the values 
obtained from the position of the maximum dissipation x" ■ We observe that J{H) exhibits barely defined crossover 
at Hi and H2 followed by more pronounced transitions at higher matching fields. It is important to note that, unlike 
the typical x'{H) measurements where several regimes are crossed while sweeping iJ, this scaling procedure allows 
one to estimate quantitatively the value of the critical current within a single dynamic regime. 

Additional information on the dynamics in this regime can be obtained by plotting x" vs. /i' = 1 + x' (Coles-Coles 
plot), as shown in figure|Slat t = 0.97 for several fields. Since in this kind of graph only dimensionless variables are 
considered, it results very suitable to analyze the data.^il5j 16] In this figure, the solid line represents the theoretical 
calculation of x"(/i') for a disk shaped sample in the Bean CS regime. For permeabilities < 0.1 a very small 
dissipation indicates the presence of a linear regime. For higher values, the similarity between the experimental 
data and the theoretical curve strongly suggest the applicability of a Bean-like model at high h. On top of that we 
observe that the maximum of x" is slightly smaller than the theoretically predicted, and it shifts towards /j,' — > 1. 
These features indicate that fiux creep is not relevant in this case, consistently with the observed weak /-dependence. 

Further confirmation that data taken at large h is well described by a Bean model can be obtained by analyzing 
the asymptotic behavior of the real part of the ac susceptibility x' at high h. Recently, Shantsev et alJ\B\ have shown 
that according to the CS model, x' cx^ h~'^ at large ft., where u depends on the J{B) dependence. For example, in 
a Bean critical state model, v = —3/2, whereas for the exponential and the Kim model, v — —3. Accordingly, we 
have found, for the whole field range studied, that v = —1.50 ± 0.15. Some of these curves are shown in the inset of 
figure El 
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FIG. 5: Experimental curves x"(m') for the data points showed in figure |1] The solid line corresponds to the theoretical 
calculation for a disk in the Bean critical state model. The inset shows the asymptotic behavior of x' for large h together with 
the linear fits from where the exponent is determined (see text). 

As we noted above, the average curvature of the pinning potential a obtained in the linear regime, together with 
the persistent current J determined in the critical state regime, provide us with complementary information about 
the principal characteristics of the pinning landscape. Usually, J is related with a through au ~ J^Oi where u is the 
pinning range. In this equation, the first term represent the restoring force due to the parabolic pinning potential 
and the second term is the Lorentz force exerted by the current J. We can use this expression in order to estimate 
the pinning range u from the two independently obtained parameters a and J. In particular, taking the values 
a{t = 0.90) « 5 X 10^ from Ref. 2], and the current density J{t = 0.90) w 5 x 10^ A/cm^ for that sample, we 
obtain a pinning range u 330 nm which is consistent with the expected range it ~ A « 400 nm. 

Finally we would like to point out that alternative models other than a critical state have been proposed for this kind 
of systems. In particular, Cooley and GrishinQ have shown that in a one-dimensional system with a periodic pinning 
array, a terraced flux proflle should be established. Similarly, later on Reichhardt et al.,'l7'| by means of molecular 
dynamic simulations, showed that in a 2D system, magnetic flux penetrate in the sample forming complex patterns 
like islands and striped domains. Although clearly, further theoretical works for higher dimensions and considering 
realistic sample geometries are needed, it is important to stress that the observed ac response can be satisfactorily 
accounted for by a simple Bean model without the necessity of invoking a multi-domain model. 
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V. CONCLUSIONS 



In summary, we have demonstrated that the Bean critical state model accurately describes the scaling observed in 
the ac susceptibility response and provides a reliable way to determine the critical current of the system as a function 
of temperature and field. The information gathered from different dynamic regimes allowed us to determine the 
averaged pinning range of the periodic landscape and estimate the boundary between the Campbell and the critical 
state regimes. 
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